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Summary: The interplay of charge, spin and lattice degrees of freedom 
is studied for quasi-one-dimensional electron and spin systems coupled to 
quantum phonons. Special emphasis is put on the influence of the lattice 
dynamics on the Peierls transition. Using exact diagonalization techniques 
the ground-state and spectral properties of the Holstein model of spinless 
fermions and of a frustrated Heisenberg model with magneto-elastic cou- 
pling are analyzed on finite chains. In the non-adiabatic regime a (T = 0) 
quantum phase transition from a gapless Luttinger- liquid/spin-fluid state 
to a gapped dimerized phase occurs at a nonzero critical value of the 
electron/spin-phonon interaction. To study the nature of the spin-Peierls 
transition at finite temperatures for the infinite system, an alternative 
Green's function approach is applied to the magnetostrictive XY model. 
With increasing phonon frequency the structure factor shows a remark- 
able crossover from soft-mode to central-peak behaviour. The results are 
discussed in relation to recent experiments on CuGeOs. 



1 Introduction 



Low dimensional electronic materials are known to be very susceptible to struc- 
tural distortions driven by the electron-phonon interaction. Probably the most 
famous one is the Peierls instability jy of one-dimensional (ID) metals: as the 
temperature is lowered the system creates a periodic variation in charge den- 
sity, called a "charge-density- wave" (CDW), by shifting the electrons and ions 
from their symmetric positions. For the half-filled band case the CDW is com- 
mensurate with the lattice and cannot slide as a whole. As a result the unit 
cell is doubled and the system has a broken-symmetry ground state. Since the 
dimerization of the lattice opens a gap at the Fermi surface the Peierls process 
transforms a metal into an insulator. Spontaneous dimerization transitions to 
a less symmetric but lower-energy configuration like those shown in Fig. [I] (left 
panel) have been found in many quasi-lD materials, such as the organic con- 
jugated polymers [e.g., (CH) Z ] and charge transfer salts [e.g., TTF(TCNQ)] or 
the inorganic blue bronzes [e.g., K0.3M0O3] and MX-chains g. 



2 



Fehske et al. 



As a generic theoretical model for such systems the ID Hubbard model is 
frequently considered, 

W c = - ^t u n la -^2t it+1 (cl a c i+la +U.c.) + U^n^riii , (1.1) 

l,(T 2,(7 i 

supplemented by a coupling to the phonon system 



according to 



SSH-type [31: tu = t ii+ i — > t(l + A?j) ft = «j - u i+ i 
Holstein-type |j: i w — > A% = t = u t , 

(1-3) 

i.e., the lattice vibrations ft interact with the electrons by modifying the electron 
hopping matrix element t y and on-site potential i ,j , respectively. In Eq. ([O ) , cJ CT 
(c ia ) creates (annihilates) a spin-cr electron at Wannier site i, and rn a = cJ la c ia . 

If there is one electron per site and the Coulomb repulsion is strong, U ^> t, 
we are in the limit of localized electrons interacting via an effective antiferromag- 
netic (AF) exchange interaction J (oc t 2 /U) and the system can be described by 
an Heisenberg Hamiltonian with magneto-elastic coupling B 

H = -jJ2( 1 + x Vi)S l S l+1 +H p , (1.4) 

i 

where = 4a^' 5 ia' witn = c il (1 - "i,-<r)- In analogy to the Peierls 

instability in ID metals, the dependence of J on the distance between neigh- 
bouring spins again gives rise to an instability, where the energy of the spin 
chain is lowered by dimerizing into an alternating pattern of weak and strong 
bonds. This so-called "spin-Peierls" (SP) transition leads to the formation of 
a singlet (dimer) ground state and there is an energy gap to elementary (mas- 
sive) spin triplet excitations being well-separated from the continuum [see Fig. |l] 
(right panel)]. Experimentally, the SP phenomenon was observed in a number 
of organic compounds, such as (TTF)(CuBDT) or MEM(TCNQ) 2 §. 

Most theoretical treatments of the Peierls instability describe the lattice de- 
grees of freedom classically. In a wide range quasi-lD metals, however, the lattice 
zero-point motion is comparable to the Peierls lattice distortion, which makes 
the rigid lattice approximation questionable 0. By any means lattice dynami- 
cal (quantum phonon) effects should be included in a theoretical analysis of the 
extraordinary transport and optical phenomena observed in Peierls-distorted 
systems ||, ||. Likewise the interest in models of spins coupled dynamically to 
phonons has increased significantly since it was recognized that the first inorganic 
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Figure 1 Schematic representation of the Peierls- and spin-Peierls scenarios. 
The left panel shows the opening of a gap A in the electronic band structure 
E(k) of an ID metal at the Fermi surface if, according to an SSH- or Holstein- 
type of electron-phonon coupling, a static lattice distortion with the new lattice 
period 2a occurs. The right panel represents the main features of the excitation 
spectrum of a Peierls-distorted ID quantum spin chain. Above the singlet (S) 
ground state at least one elementary excitation, corresponding e.g. to a triplet 
(T) bound state, is split from the continuum (C). 



SP compound CuGeC>3 |tC]] shows no clear separation between the magnetic 
and phononic energy scales: the two (weakly dispersive optical) Peierls-active 
phonon modes have frequencies wo,i — J and wo,2 — 2 J |ll], [12]. Moreover, 
in contrast to the organic SP materials, no phonon-softening is observed at the 
SP transition The SP physics in CuGeC>3 is therefore in the non-adiabatic 
regime. 

Motivated by this situation, in this report, we study the perhaps minimal 
microscopic models capable of describing the Peierls and spin-Peierls transi- 
tions in ID systems by the use of numerical techniques allowing an essentially 
exact treatment of both subsystems, electrons/spins and phonons, at a fully 
quantum-mechanically level. In the weak-coupling regime, the random-phase- 
approximation (RPA) approach is shown to be consistent not only with phonon 
softening but also with phonon hardening at the SP transition, as observed, e.g., 
in (TTF)(CuBDT) and CuGe0 3 , respectively. 
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2 Luttinger-liquid vs. charge-density-wave behaviour 



First we consider the ID Hubbard model, Eq. (1.1), at quarter filling and con- 
fine ourselves to the case of spinless fermions for simplicity. This model is of 
physical relevance in the strong interaction limit U — ► oo and is particularly 
interesting because a quantum phase transition from a Luttinger liquid (LL) to 
a CDW phase occurs at a finite electron-phonon interaction, demonstrating that 
the quantum phonon fluctuations destroy the dimerized ground state for weak 
electron-phonon couplings [jli], p"5| . 

Rescaling H — > H/t = H c + H p + H cp and setting 



A = g^2Kuj with uj% = K/M (2.5) 
the Holstcin Hamiltonian in a particle-hole symmetric notation is given by 

We = -$>i c «+i+4+icJ, «p = w oE(^, + I). (2-6) 



and 



H e - P = -goj £>t + K) K - \) ■ ( 2 - 7 ) 



In Eqs. (2.6) and (2.7), bp annihilates (creates) a dispersionsless Einstein phonon 
of frequency loq coupled to the local electron density rii = cjcj (K = 1, and all 
energies are measured in units of t) . Note that for the Holstein model the dimen- 
sionless electron-phonon coupling constant g — y/ £ p /co>o is directly related to the 
familiar polaron shift e p being a second natural measure of the strengths of the 
electron-phonon interaction. Both parameters are necessary in order to charac- 
terize the weak (e p <€. 1) and strong coupling {e p 3> 1 and g » 1) situations in 
the adiabatic (loq <C 1) and anti-adiabatic regimes (loq 3> 1), respectively. 

Previous results for the ground-state phase diagram of the Holstein model 
at half-filling obtained by WL QMC |4) and GF MC @ simulations showed 
significant discrepancies in the region of small uiq (0 < uiq Only very recently 
Bursill et al. jlBj provided more reliable information from level crossings in their 
DMRG data. Applying a new optimized phonon approach for the diagonalization 
of coupled electron/spin-phonon systems p7| , we consider the Holstein model on 
chains of even length with up to 10 sites and periodic (antipcriodic) boundary 
conditions if there is an odd (even) number of fermions in the system. The 
resulting phase diagram is shown in Fig. || over a wide range of frequencies 
and coupling strengths. For small g the system is a metal, more precisely a 
Luttinger liquid with parameters that vary with the coupling (see below). For 
large g the system has an energy gap and develops true long-range CDW order 
in the thermodynamic limit. The phase boundary obtained with our optimized 
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Figure 2 Ground-state phase diagram of the ID Holstein model of spinless 
fermions at half- filling (N e = N/2), showing the boundary between the Luttinger 
liquid (LL) and charge-density-wave (CDW) states obtained by exact diagonal- 
ization (ED) and density matrix renormalization group (DMRG) mM approaches. 
The dashed line gives the asymptotic result for the XXZ model. Lett insets show 
the LL parameters u p and K p as a function of the electron-phonon g in the 
metallic regime; right insets display for a six-site chain the regular part of the 
optical conductivity a re9 (ui) (dotted lines) and the integrated spectral weight 
S re3 (uj) = J " dLu'o re9 {uj') (solid lines) in the CDW region. 



phonon diagonalization methods agrees with recent DMRG results jl5|. In the 
adiabatic limit luq = 0, the critical coupling converges to zero, as expected for the 
Holstein Hamiltonian, Eqs. ( |l.l| )-(1.3), with M — > oo. For < ujo ^ 1 the accurate 
determination of g c is somewhat difficult. On the other hand, in the strong- 
coupling anti-adiabatic regime, the half-filled Holstein model can be transformed 
to the exactly soluble XXZ (small polaron) model jl4|, |l8) 



^XXZ = N {2ujo g2 ^ y ^ c _ g , £ [ {S + S - +i + g-s+j _ y 2 ^gig 



4+1 

(2.8) 
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using second order perturbation theory with respect to t. Here V%(jg = 
2e~ 2g uJq 1 J2 s ^ Q (2g 2 ) s /(ssl), and the (Kosterlitz-Thouless) phase transition line 

is given by the condition V^a, g 2 )e g2 /2 — 1 (long-dashed curve in Fig. ||). 

Let us now characterize the LL and CDW phases in some more detail. Accord- 
ing to Haldane's Luttinger liquid conjecture ID gapless systems of inter- 
acting fermions should belong to the same universality class as the Tomonaga- 
Luttinger model. As stated above, the Holstein system is gapless for small enough 
coupling g. Thus it is obvious to prove, following the lines of approach to the 
problem by McKenzie et al. jig ], whether our Lanczos data shows a finite-size 
scaling like a Luttinger liquid. For a LL of spinlcss fermions, the ground-state 
energy Eq (N) of a finite system of N sites scales to leading order as 

N °° 67V 2 ' [ ' 

where £oo denotes the ground-state energy per site for the infinite system and u p 
is the velocity of the charge excitations. If E±\(N) is the ground-state energy 
with ±1 fermions away from half filling, to leading order the scaling should be 

77-7 / 

E ±1 {N)-E {N) = W ^. (2.10) 

K p is the renormalized effective coupling (stiffness) constant. The left insets of 
Fig. H show the LL parameters as a function of g in the region, where the scaling 



relations Eqs. ( |2.9| ) an d ( 2.1C ) hold. A very interesting result is the changing 
character of the interaction below u>q ~ 1. For small frequencies the effective 
fcrmion-fcrmion interaction is attractive, while it is repulsive for large frequen- 
cies, where the system forms a polaronic metal with strongly reduced kinetic 
energy JT^ ] . Self-evidently there is a transition line in between, where the model 
describes "free" particles in lowest order. 

In the CDW state extremely valuable information about the low-energy exci- 
tations can be obtained from the behaviour of the optical conductivity. The real 
part of <j{uj) contains two contributions, the (coherent) Drude part at uj = and 
a so-called "regular term", a re9 (u), due to finite- frequency dissipative optical 
transitions to excited quasiparticle states. In spectral representation (T = 0), 
the regular part takes the form 

gre>) = E 5[u _ {Em _ Eo)] t (2 . n) 

m>0 m 

where a re9 (uj) is given in units of ire 2 and we have omitted an 1/iV prefac- 
tor. The evaluation of dynamical correlation functions, such as Eq. ( |2.11| ), can 
be carried out by means of very efficient and numerically stable Chebyshev re- 
cursion and maximum entropy algorithms EG]. Clearly the optical absorption 
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spectrum in the strong EP coupling regime is quite different from that in the 
LL phase (cf. Ref. |^D). ft can be interpreted in terms of strong electron-phonon 
correlations and corroborates the CDW picture. Since for g > g c the electronic 
band structure is gapped we expect that the low-energy gap feature survives in 
the thermodynamic limit. In the adiabatic region (upper right inset), the broad 
optical absorption band is produced by a single-particle excitation accompanied 
by multi-phonon absorptions and is basically related to the lowest unoccupied 
state of the upper band of the CDW insulator. The lineshape of a re9 (uj) reflects 
the phonon distribution in the ground state. The most striking feature is the 
large spectral weight contained in the incoherent part of optical conductivity. 
Moreover, employing the f-sum rule for the optical conductivity p^] and taking 
into account the behaviour of the kinetic energy (oc u p ) as function of g, we 
found that in the metallic LL and insulating CDW phases nearly all the spec- 
tral weight is contained in the coherent (Drude) and incoherent (regular) part 
of Re <j{uS), respectively. As stated above, in the anti-adiabatic regime the LL 
phase is basically a polaronic metal, i.e., the electrons will be heavily dressed 
by phonons. Since the renormalizcd coherent bandwidth of the polaron band is 
extremely small, the hnite-size gaps in the band structure are reduced as well. 
Therefore, the gap occurring in the CDW state (Aqdw ~ %>) may be identified 
with the optical absorption threshold (see lower right inset). 



3 Non-adiabatic approach to the spin-Peierls transition 
3.1 Exact diagonalization results for T=0 

In spite of the experimental fact that a realistic modeling of the inorganic SP 
compound CuGeC>3 should include the phonon dynamics, previous theoretical 
studies have commonly adopted an alternating and frustrated AF Heisenberg 
spin chain model p3| 

^static = [(1 + Si-l^Sji+i + a Sji+2] (3.12) 

i 

with a static dimerization parameter 5, thus representing the extreme adiabatic 
limit of a SP chain (in this section all energies are given in units of J), a deter- 
mines the strength of the frustrating AF next-nearest-neighbour coupling. The 



spin model (3.12) contains two independent mechanisms for spin gap formation. 
At S = and for a < a c the ground state is a spin liquid and the elementary 
excitations are massless spinons [Q. The critical value of frustration a c — 0.241 
was accurately determined by numerical studies |2^, |25| . For a > a c the ground 
state is spontaneously dimerized, the spectrum acquires a gap, and the elemen- 



tary excitations are massive spinons [E6, 27 1. On the other hand for any finite S 
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the singlet ground state of the model ( 3.12 ) is also dimerized, but the elementary 



excitation is a massive magnon 



24 



_28| . A comprehensive study of the spectral 
properties of the model ( 3.1 2| ) in terms of the spin dynamical structure factor 
has been carried out by Yokoyama and Saiga |29| . 

From the magnetic properties of the uniform phase J ~ 160 K and a = 0.36 
have been estimated for CuGeC>3 [[30). However, if one attempts to re produce 
the observed spin gap A ST ~ 2.1 mcV within the static model ( 3.12 ), a very 
small value of S ~ 1.2% results. From the uniaxial pressure derivatives of the 
exchange coupling J [3l| ] and the structural distortion in the dimerized phase ]ll[ | 
a minimum magnetic dimerization of about 4% is obtained, incompatible with 
an adiabatic approach to the SP transition. 

The simplest model that maintains the full quantum dynamics of the lattice 
vibrations may be obtained from Eq. ( 3.1 2| ) by replacing (— l) l S — > gujo(b\ + b i ); 



i+1 



■aSiSi +2 ), Hp = uj q b\bj , 



(3.13) 



Hl p =9^^2(4 + ^)^+1. (3.14) 

i 

Recently it was shown that such a dynamical spin-phonon model describes the 
general features of the magnetic excitation spectrum of CuGeOa |3^, |3^, |34[ ]. 
Here we focus on the behaviour of the lattice dimerization which can be found 
from the displacement structure factor at wave number q = tt 

5 2 = ^yj2 C ^ ARi ~ Ri) with C?? = ((b i + bt)(b j+ b})). (3.15) 

The alternating structure of the correlation function Cu as shown in the inset 
Fig. H (a) implies the Peierls formation of short and long bonds and thus al- 
ternating strong and weak AF exchange interactions, i.e., a dimerized ground 
state. The structure is enhanced (weakened) increasing the spin-phonon coupling 
(phonon frequency) |S2| . As in the ordinary Peierls phenomenon, a finite dimer- 
ization S > necessarily leads to a gap A ST in the magnetic excitation spectrum. 
For evaluating the relation between the dimerization and the resulting magni- 
tude of the spin triplet excitation gap, we keep the phonon frequency fixed, vary 
the coupling strength g and calculate for each parameter set A ST = Eq — Eq and 
the dimerization 6 from Eq. ( |3.15| ). ED results for A ST obtained for the static 
spin-only model ( |3.13| ) and the quantum phonon model (3.14) are compared in 



the main part of Fig. |3| (a). For vanishing dimerization, i.e. in the absence of 
any spin-phonon coupling (<? = 0) the results for the static and the dynamic 
model naturally agree (note that for 5 = there remains a spin excitation gap 



due to the frustration driven singlet dimer ordering). The dynamic model (3.14) 



partially resolves the A ST — S conflict we are faced within the static approach, 
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Figure 3 Spin gap (a) and dimerization (b) in the frustrated Heisenberg spin 
chain with dynamical spin-phonon coupling. 



because the dimerization S grows with the phonon frequency at fixed A . Thus 
matching A ST to the neutron scattering data, a larger S ~ 5% may result. 

Figure |^ (b) shows the dimerization as a function of the spin-phonon interac- 
tion strength. As for the Holstein model the precise determination of the critical 
coupling g c is difficult for phonons in the adiabatic regime u>q < 0.5. On the 
contrary, in the non-adiabatic regime the SP transition takes place at g c ~ 1 
nearly irrespective of ojq. Most remarkably, below (above) g c the dimerization 
decreases (increases) with increasing lattice size N (see inset), indicating that 
the infinite system exhibits a true phase transition. 



3.2 Anti-adiabatic limit: mapping to an effective magnetic problem 

Figure ^ (b) implies that more insight into the dynamic spin-phonon model 
can be obtained by considering the limit of large phonon frequencies. One can 
then integrate out the lattice degrees of freedom in order to derive an effective 
spin Hamiltonian which cover the dynamical effect of phonons in a approximate 
way. Technically this can be achieved by a variety of methods, such as standard 
perturbation theory f35| , continuous (flow-equation based) [ |36| or variational 
unitary transformations fl37f. 



In what follows we consider besides the magneto-elastic interaction, Eq. ( 3.14 ) 



another type of spin-phonon coupling |35 



Hll v = <?o> £>l + b i)(SiS i+1 - SiSi-x) , (3.16) 



where the AF exchange integral varies linearly with the difference between the 
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phonon amplitudes on neighbouring sites, and perform a Schricffcr- Wolff trans- 
formation |39), TL = exp(iS)H exp(— S), with 

5 I = fl ^(i!-t 1 )S4 1 and S n = gJ2(bt-b i )(S i § i+1 -S i S i _ 1 ), (3.17) 

i i 

respectively. In contrast to electron-phonon systems with Holstein coupling, 
where the (Lang-Firsov) transformation similar to exp(S') completely removes 
the electron-phonon interactio n ter m (7i c _ p ), applying the unitary transforma- 
tion exp(S) to H with Eqs. (fig) and |H|), we obtain an infinite series o:' 



terms, which can not be summed up to a simple expression. To derive an effec- 
tive spin model we now take the average over the (transformed) phonon vacuum. 
The resulting spin Hamiltonian TL c s, s — (0 p |7i|0 p } contains longer than next- 
nearest-neighbour ranged Heisenberg interactions as well as numerous multi-spin 
couplings. To a good approximation we can neglect them and obtain (cf. J37J) 

^cft - + a c ^SA +2 ) (3.18) 



with 



1 a + ff 2 (l - 2a)/2 + 3/^/16 - 37g 4 (l - 2a)/96 

a ° ff l + ff 2 ^o/2-3 2 (l-a)/2-3 5 4 w /8 + 3 4 (28 - 37a)/96' l ' J 



+ g 2 u Q /2 + g 2 (3 - 5a)/ 2 + 3 5 4 w /48 - 3 4 (75 - 124a)/24 
1 + .g 2 w - 3g 2 (l - a)/2 - 9g 4 LU /8 + ,g 4 (59 - 75a)/24 



(3.20) 



That is the integration over the phonon subsystem yields an additional frustrat- 
ing next-nearest-neighbour exchange interaction. Therefore, without any explicit 
frustration a, the effective frustration a e ff due to the phonons can lead to a gap 
in the energy spectrum and to spontaneous dimerization. This effect is most 
important in the anti-adiabatic frequency range |3jj ] . 

In order to determine the critical line in the the a — g plane indicating the 
transition from the gapless to the gapped ground state, we use the level crossing 
criterion j2^, |23|, [37| |38| for the lowest singlet and triplet excitations which 
become degenerate at a c (N). At a c the finite-size corrections a c (N) — a c (oo) ~ 
7V~ 2 are small. The resulting phase diagrams are displayed in Fig. [|. We find 
that the phase boundary of the original quantum spin-phonon and effective spin 



models agree surprisingly well, where for the effective model (3. IE) a c (g) can 
be obtained with high accuracy on local workstations if N < 20. For spin- 
phonon coupling of type I the critical curve exhibits a remarkable upturn before 
crossing the abscissa; i.e., the frustration is suppressed for small spin-phonon 
coupling, but over-critical for strong coupling. It is this feature which makes it 
necessary to expand TL up to fourth order in g to approximate H in a correct 
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Figure 4 Phase d iagra m of the original spin-phonon model (3.13) and the ef- 
fective spin model (3.16) with spin-phonon couplings of types I (a) and II (b). 



way. A second order theory is not capable to describe the observed critical line. 
The upturn of a c (g) at small g was reproduced quite recently by linked series 
expansion techniques pof . By contrast a c (g) is a monotonous decreasing function 
of g for the coupling of type II. It appears that one would get the same shape 
also for a second order theory. However, to enlarge the application area of our 
approximation taking into account higher order contributions is still appropriate. 
For a = 0, differently from the coupling case I, g c tends to zero in the anti- 
adiabatic limit loq — > oc. While the q = and the q = ir phonon mode compete 
in the case of coupling I, allowing for a stable gapless phase up to a critical 
<?, there is no interaction with the q = mode in flS. Therefore the q = it 
mode induces long ranged exchange more efficiently, leading to a vanishing g c 
for uq — » oo. 



3.3 RPA approach at T ^ 0: Soft-mode vs. central-peak behaviour 

The absence of a soft phonon mode at the displacive SP in CuGeC>3 has been 
puzzling for a long time, since the behaviour was believed to be inconsistent 
with the standard Cross-Fisher RPA approach to SP transitions JO]. Even 
worse, the relevant phonon modes harden by about 5% with decreasing temper- 
ature |Q, pointing to a central peak scenario. A way to reconcile these results 
within the framework of the Cross-Fisher theory has been proposed by Gros and 
Werner [43)] , with the result that a complete phonon softening occurs only for 
bare phonon frequencies less than a critical value (luq < wo, c = 2.2T C , where T c 
is the SP transition temperature. For higher luq a central peak develops reaching 
T c from above. For the magnetostrictive Heiscnberg spin chain model, however, 
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it was not possible to analyze the complete pole structure of phonon spectral 
function. 

To test this scenario we restrict ourselves to the simpler XY model 

nf Y = J2(SfSf +1 + sys^ +1 ), h p = J2vI + ^- u ^ 2 ( 3 - 21 ) 



with the spin-phonon coupling term 



W s -p = \ - u t+ i)(S+Sr +1 + S r S + +1 ) . (3.22) 

i 

The above magnetostrictive XY model has been proposed as the minimal model 
to describe the SP transition [j44| . Performing a Jordan- Wigner transforma- 
tion |45f it can be mapped onto a model of spinless fermions only interacting 
with the phonon system H p = ^2 q ^J q b q b q : 

Hf Y - = -^cos(fc)4c fc , (3.23) 

k 

H s _p - Hi w p = gu; J2v(k,qKb q + bl q )c{ Ck _ q , (3.24) 

fe . g 

where 

i 

V(k, q) = -i 2 [sin(fc -q)- sin(fc)] . (3.25) 

with w g = 2u)$ sin(g/2), i.e = 2wo- 

Determining the equation of motion for the Matsubara phonon Green's func- 
tion D(q,ioj n ) within the scheme worked out by Bennett and Pytte for the 
magnetostrictive Heisenberg model jl6| , the main advantage is that the higher- 
order spin-spin correlation functions can be calculated without further approx- 
imations for the XY-case. The retarded Green's function D ret (q, ui) is obtained 
from D(q, iuj n ) by analytical continuation ito n —> lu. In the uniform phase above 
T c , the RPA propagator of the q = 7r-phonon, being responsible for the phase 
transition, then results as 



D le \n,uj) = - 2 AUJ0 -, (3.26) 

UT — 4uJq — 2ldoP(TT, UJ) 



r dk sin 2 (fc)tanh[f cos(fc)] 

P(w,u>) = -(2guj Q ) / — — . 3.27) 

7T cj + 2cos(fc) 
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Figure 5 Dynamical structure factor for the magnetostrictive XY model. 

A structural instability occurs if D Tet (q, ui) exhibits a pole at a certain q- value for 
uj = 0. In that case an excitation with arbitrary low energy is possible and the 
lattice becomes unstable with respect to a (static) deformation with wave num- 
ber q. From this condition the equation for the inverse transition temperature 
(3 C = 1/T C easily follows as 

r sin 2 (fc)tanh[%cos(fc)l 

1 = k dk U i-f Ui, 3.28 

Jo cos(fc) 

where k = g 2 ujo/Tr. A more detailed study of the complete pole structure of 
D Tet (q, uj) in both the undimerized and dimerized phases is presented in Ref. j47|, 
where also the ultrasound properties of the system were discussed. Here we focus 
on the temperature dependence of the dynamical structure factor above T c , 

C M = -inra I ^ et ^ +iS \ (3.29) 

which is displayed in Fig. |5| for two characteristic phonon frequencies ujq = 
0.25 and ui = 1.05, corresponding to the soft-mode (a) and central-peak (b) 
regimes, respectively. For small bare phonon frequencies the maximum in C(n, ui) 
is located around lo = lu^ at high temperatures and moves with decreasing 
temperatures to lower frequencies until a real divergence of C(it,uj) appears 
for T = T c at uj = 0. For large bare phonon frequencies a completely different 
behaviour is found. Now the high-temperature phonon peak stays around to = uj^ 
and even hardens in some degree as T — > T c . At the same time a peak structure 
evolves at q = which becomes divergent at T = T c . This so-called central peak 
corresponds to a new collective magneto-elastic excitation of the coupled spin- 
phonon system, which occurs at the SP transition and moves to higher energies 



as the temperature is lowered further 47 
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4 Conclusions 

In this report the lattice dynamical effects on the Peierls transition in coupled 
electron/spin-phonon systems were discussed. Our approach was based on several 
generic model Hamiltonians obtained from the one-dimensional Holstcin/SSH 
Hubbard model in important limiting cases. Applying both basically exact nu- 
merical and approximative analytical methods we are able to calculate both 
ground-state and spectral properties of these simplified models for the complete 
range of model parameters. The focus, however, was on non-adiabatic effects due 
to the phonon dynamics. The presented results confirm previous findings that 
quantum phonon fluctuations destroy the Peierls distorted state at sufficiently 
weak electron/spin-phonon interactions. 

For the spinless fermion model this means that for weak electron-phonon cou- 
plings the system resides in a metallic (gapless) phase described by two non- 
universal Luttingcr liquid parameters. The renormalized charge velocity and the 
correlation exponent are obtained from finite-size scaling relations, fulfilled with 
great accuracy. The Luttinger liquid phase splits in an attractive and repul- 
sive regime at low and high phonon frequencies, respectively. Here the polaronic 
metal, realized for repulsive interactions, is characterized by a strongly reduced 
mobility of the charge carriers. Increasing the electron-phonon coupling, the 
crossover between Luttinger liquid and charge-density-wave behaviour is found 
in excellent agreement with very recent DMRG results. The transition to the 
CDW state is accompanied by significant changes in the optical response of the 
system. Most notably seems to be the substantial spectral weight transfer from 
the Drude to the regular (incoherent) part of the optical conductivity, indicating 
the increasing importance of inelastic scattering processes in the CDW (Peierls 
distorted) regime. 

As a simple model for a spin-Peierls system we considered a (frustrated) 
Hciscnberg spin chain model supplemented by a coupling to optical phonons with 
frequencies comparable to the magnetic exchange coupling, which is, e.g.. the 
relevant limit for the spin-Peierls compound CuGeOs. The magnetic excitations 
inherently include a local lattice distortion requiring a multi-phonon-mode treat- 
ment of the lattice degrees of freedom. When compared to the static model of an 
alternating, dimerized spin chain the magnetic properties are strongly renormal- 
ized due to the coupled spin and lattice dynamics. For the quantum spin-phonon 
model the dimerization dependence of the spin triplet excitation gap is found 
to be in qualitative agreement with experiment. In the anti-adiabatic phonon 
frequency range we used the concept of unitary transformations in order to inte- 
grate out phonon degrees of freedom, where the spin-phonon interaction results 
in an additional frustration in the effective spin-only model. For two types of spin- 
phonon couplings, the critical lines separating gapless from gapped ground states 
are found to agree for the original spin-phonon and effective spin models. Finally 
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the reliability of the RPA approach to SP transitions was assessed in terms of an 
XY model with additional magneto-elastic interaction. With increasing phonon 
frequency the dynamical displacement structure factor shows a crossover from 
a soft-mode to central-peak behaviour, normally linked to displacive and order- 
disorder types of transition respectively. The central peak predicted to appear 
at T c for large bare phonon frequency corresponds not to a real phononic but 
rather to a new magneto-elastic excitation. 
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